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ABSTRACT 

C$ ■ We use the standard, adiabatic shell evolution to predict the differential size distribu- 

tion N(R) for populations of OB superbubbles in a uniform ISM. Assuming that shell 
growth stalls upon pressure equilibrium with the ambient ISM, we derive N(R) for sim- 
l> I pie cases of superbubble creation rate and mechanical luminosity function (MLF) . For 

constant creation and an MLF 0(L) oc L' 13 , we find that N(R) oc i? 1_2/3 for R < R e , 
and N(R) oc R 4 ~ 5 * for R > R e , where the characteristic radius R e ~ 1300 pc for 
J> ■ typical ISM parameters. For R < R e , N(R) is dominated by stalled objects, while for 

R > R e it is dominated by growing objects. The relation N(R) oc R 1 ^ 213 appears to be 
quite robust, and also results from momentum-conserving shell evolution. We predict 
\ a peak in N(R) corresponding to individual SNRs, and suggest that the contribution 

■ of Type la SNRs should be apparent in the observed form of N(R). We present ex- 
| pressions for the porosity parameters, Q2D and Q3D, derived from our analysis. Q2D 

Q\ ■ is dominated by the largest superbubbles for (3 < 2 and individual SNRs for (3 > 2, 

whereas Q3D is normally dominated by the few largest shells. 
O We examine evolutionary effects on the HII region luminosity function (HII LF), 

in order to estimate (5. We find that for a nebular luminosity fading with time t, C oc 
; 1 ■ t~ v , there is a minimum observed slope a m ; n for the HII LFs. Empirical measurements 

I all show a > a m m> therefore implying that usually we may take (3 — a. We also find 

d • that if nebular luminosity is instantaneously extinguished at some given age, rather 

than continuously fading, no a m i n will be observed. 

Comparison with the largely complete HI hole catalog for the SMC shows sur- 
prising agreement in the predicted and observed slope of N(R). This suggests that 

■ no other fundamental process is needed to explain the size distribution of shells in 

the SMC. Further comparison with largely incomplete HI data for M31, M33, and 
Holmberg II also shows agreement in the slopes, but perhaps hinting at systematic 
differences between spiral and Im galaxies. We estimate porosities that are substan- 
tially < 1 for all of the galaxies except Holmberg II, for which we obtain values il 1. 
Most of these galaxies therefore may not be strongly dominated by a hot interstellar 
component. However, porosity results for the Galaxy remain inconclusive with the 
available data. 

Key words: ISM: bubbles — ISM: general — HII regions — ISM: structure — 
supernova remnants — galaxies: ISM — galaxies: individual: SMC 



1 INTRODUCTION 

The evolution of superbubbles created by the stellar winds and supernovae of OB associations is one of the primary processes 
that determines the structure and energetics of the interstellar medium (ISM). All components of the diffuse interstellar 
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medium, including warm and cold HI as well as the warm ionized medium (WIM), are thought to have large fractions of 
their volume consisting largely of superbubble walls {e.g., Heiles 1984; Kennicutt etal. 1995). The restructuring of cool gas 
will in turn influence the collapse of molecular clouds and attendant star formation {e.g., McCray & Kafatos 1987). Likewise, 
the hot ionized medium (HIM) is believed to originate within these superbubbles and Type la supernova remnants (SNRs), 
hence the evolution of these structures determines whether and how this gas is released into the general ISM. Models for the 
ISM are strongly dependent on whether the coronal gas is the pervasive, dominant component, as envisioned by e.g., McKee 
& Ostriker (1977), or whether it plays a lesser role as favored by e.g., Slavin & Cox (1993), returning to a paradigm closer to 
the two-phase ISM of Field, Goldsmith, & Habing (1969). 

The standard model for understanding the evolution of superbubbles and other shell structures is that of an adiabatic, 
pressure-driven bubble, with continuous wind energy injection {e.g., Weaver etal. 1977; Dyson 1977; Pikel'ner 1968). It has 
generally been assumed that in later stages, successive supernovae (SNe) can be treated as an approximate continuous energy 
injection, and are able to power the growth of supergiant shells {e.g., Mac Low & McCray 1988). Such supergiant shells, 
having radii of hundreds of pc, could then blow out of the galactic disk, releasing hot gas into the halo, possibly in a galactic 
fountain cycle. 

However, recent evidence suggests that superbubble evolution is often not as simple as described by the standard model. 
Many superbubbles exhibit expansion velocities that are too high to be consistent with the standard evolution, and a correlation 
of bright X-ray emission with many of these objects suggests acceleration by SNR impacts (Oey 1996) . In such cases, not all 
the available SN energy will be thermalized to power the shell expansion, thus the detailed effect that discrete SNe have on the 
long-term shell evolution remains unclear. Likewise, there is evidence that the growth rate of the shells may be overestimated, 
perhaps owing to an overestimate in input wind power predicted by assumed stellar mass-loss rates (Oey 1996; Garcfa-Segura 
& Mac Low 1995; Drissen etal. 1995). Therefore much remains unclear about the details of superbubble evolution. 

A demonstrative example is the superbubble DEM 152 (N44) in the Large Magellanic Cloud. ROSAT observations 
of this object clearly show hot gas escaping through a "blowout" feature (Magnier etal. 1996), and this object also exhibits 
kinematics discrepant from the standard model (Oey & Massey 1995). If the ISM has low-density channels that encroach upon 
the superbubble environment, such evolutionary disruptions could signal a significant release of hot gas into the surroundings, 
and not necessarily out of the plane of galactic disks. The details of superbubble evolution thus have a direct consequence for 
the relative volumes of coronal gas in galactic disks and halos. Examples such as DEM 152 demonstrate the uncertainties in 
our understanding of superbubble evolution, and the importance to the global ISM. 

One approach to test the assumed evolution is an examination of the superbubble size distribution in galaxies. This can 
be predicted from the evolutionary model in combination with a given production rate and mechanical luminosity function. A 
subsequent comparison with observed size distributions can then yield insight on the evolutionary assumptions and possible 
consequences for the inferred structure of the ISM. The cumulative size distribution of radiative SNRs has long been used 
as a diagnostic of SN parameters {e.g., Hughes, Helfand, & Kahn 1984), but it can also be used to test evolutionary models 
for SNRs and superbubbles. An earlier investigation with a foray into this technique was carried out for SNRs by Cioffi & 
Shull (1991), who pointed out that the differential size distribution is a much more powerful diagnostic than the cumulative 
size distribution. In this work, we present a prediction for the differential size distribution for superbubbles and compare the 
results with HI observations of nearby galaxies. We will then examine the results with regard to the structure and porosity 
of the ISM. 

We will consider a rather rudimentary model for superbubble evolution, with the aim of identifying dominant effects, 
rather than precisely reproducing observations. For example, the size distribution of Sedov- Taylor SNRs is predicted to increase 
as N{R) oc R 5 / 2 (Hughes et at 1984), for uniform ISM and stellar properties, whereas H I hole distributions in galaxies all show 
size distributions that decrease with shell radius R. Our objective is to derive gross general properties that would be expected 
as a result of the simplest practical application of the standard model. Agreements and disagreements with observations 
will then be instructive in understanding the features of the model that dominate the resultant size distribution and shell 
evolution, with attendant consequences for the structure of the ISM. 

The paper is organized as follows. In §2, we work out general expressions for the assumed shell evolution and derived size 
distributions, for three combinations of shell creation and mechanical power input: (a) continuous creation, single luminosity; 
(b) single burst, luminosity spectrum; and (c) continuous creation, luminosity spectrum. In §3, we evaluate each of the solutions 
for the standard, adiabatic shell evolution. We also briefly examine the standard, momentum-conserving evolution. To obtain 
the slope of the mechanical luminosity spectrum, we examine evolutionary effects in the HII region luminosity function in §4. 
In §5, we compare our predictions with HI observations in four nearby galaxies, with attention to the slope, normalization, 
and peak of the size distributions. We also estimate supernova rates based on our analysis. The consequences are discussed in 
§6, where we examine the regimes of applicability for our predictions. Here we discuss the role of processes likely to modify 
the size distributions, and the role of alternate shell creation mechanisms. These effects are likely to be dependent on galaxy 
type. We also derive expressions for the ISM porosities resulting from our analysis, and compute porosities for the four nearby 
galaxies and the Milky Way. Our results are summarized in §7. 
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2 GENERAL ANALYTIC EXPRESSIONS 
2.1 Assumptions 

Superbubble expansion is presumed to be powered by the stellar winds, and especially, SNe of the parent OB association. 
It has been customary to represent the SN power as a continuous energy injection analogous to a wind (e.g., Mac Low & 
McCray 1988), in which over half the SN energy is thermalized to drive the shell expansion. We will adopt this standard, 
adiabatic representation (e.g., Weaver etal. 1977; Dyson 1977). For the purpose of comparing with HI "hole" data and ISM 
porosities, we are particularly interested in the superbubble cavity, which in this thin-shell approximation is considered to 
evolve identically to the outer radius. We also assume that the mechanical luminosity L(t), dominated by the SNe, remains 
constant as a function of time t. As shown by Shull & Saken (1995), this is a reasonable approximation for realistic slopes of 
the stellar initial mass function (IMF). This constant L(t) is then assumed to continue until the lowest-mass SN progenitors 
expire at t = t e , which corresponds to a stellar mass of about 8Mq, or a period of roughly t e = 40 Myr. A variation in input 
power may occur in the initial stages (Ss 3 — 4 Myr) while L is dominated by the stellar wind phase (e.g., Shull & Saken 
1995), but this stage is short relative to t e , the total period of energy input. Leitherer & Heckman (1995) and Ferriere (1995) 
show that a starburst over this period does not substantially alter this approximation. A uniform ambient medium is assumed 
throughout, although in §6 we will discuss effects expected from disk galaxy gas distributions. 

As will become apparent below, the endstage of shell evolution is vitally important to this analysis, but is fraught 
with uncertainties. We will consider the following simple model, keeping in mind our motivation of testing the simplest- 
case standard evolution. In most cases, eventually the superbubbles' internal pressure Pi Pq, the ambient pressure of the 
interstellar medium, while still at ages t < t e . We will consider that for such objects, the superbubble growth stalls when 
Pi — Po- At this stage, radiative losses are thought to become important, confining the shell growth. Simulations by Garcia- 
Segura & Franco (1996), for example, show that the growth of the cavity radius does stall at times close to the sonic point. 
We will assume that the superbubbles survive at the stall radius Rf until the input power stops at time t e . In principle, the 
shells will begin to disintegrate, owing to ambient random motions. However, the disintegration timescale may be fairly long; 
for example, it is Rf/v rnls =10 Myr for ambient rms velocity v lrns — 5 km s and Rf = 50 pc. Furthermore, the process of 
disintegration will still leave detectable cavities, so it is unclear how this effect can be quantified. We note that if the breakup 
of objects into smaller subunits occurs such that the ratio of subunit sizes is universal for all objects, then a power law size 
distribution remains unaffected (Clarke 1996). Our analysis will therefore not consider elimination of shells at ages t < t e , 
nor differential elimination among the shell population. Finally, an important feature regarding stalling is the existence of a 
monotonic correspondence between L and Rf, for a uniform ISM. 

Objects that never achieve pressure equilibrium with the ISM are simply assumed to grow until t e . Finally, all objects 
are assumed to survive for an additional nominal period t s <Ct, which is the same for all shells. 

Within individual clusters, we take the star formation to occur in a single, instantaneous burst. Massey etal. (1995a, b) 
find that in general the duration of star formation is ^ 3 Myr for associations in the Magellanic Clouds and the Galaxy, which 
again is short relative to t e . The spectrum of associated mechanical luminosities is a critical parameter; we will consider both 
a single, global value for L, and a power- law luminosity function (LF): 

^ = <f>(L)=AL-? ! , (1) 

normalized such that J 4>(L) dL — 1. In §4 we will discuss the form of </>(!/), and in particular, the value of (5; our analysis 
of 4>(L) will assume the stellar IMF to be constant. We will also consider a constant cluster formation rate ip(i), and an 
instantaneous burst of cluster formation. 

We first present the general expressions for the superbubble size distribution in the case that the shells grow according 
to the generic law: 



R — min 



(R(L,t), Rf(L)^j . (2) 



where Rf(L) is the radius at which a superbubble stalls, which is dependent only on the input L, for uniform ISM parameters. 
We now consider the production of superbubbles for the simple cases of luminosity distribution 4>(L) an d formation rate tp(t) 
mentioned above, with respect to the stalling and continuous-growth evolutionary schemes. In each case, we derive N(R), 
defined such that the number of superbubbles with radii in the range R to R + dR is N(R) dR. 



2.2 Continuous Creation, Single Luminosity 

If superbubbles are generated at a constant rate ip, then the number of growing shells with radii in the range RtoR + dR is 
equal to ip dt where dt is the time interval for creation corresponding to this radial range. Thus, the differential size distribution 
is, 
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iWfl) = vf^J ■ (3) 

[Note that partial derivatives with respect to t and L are evaluated at constant L and t respectively.] 

Recalling the correspondence between L and stall radius Rf, the single- valued cj>(L) — Lo will therefore yield a single- 
valued iVstaii = N(Rf) for stalled objects: 

N stall (R) = ^(t e -tf(L )^ -S^R-RfiLo)^ , (4) 
where i/(Z/o) is the age at which a shell powered by Lo stalls. 



2.3 Single Burst, Luminosity Spectrum 

If N), superbubbles are created in an instantaneous burst, then after time t, the number of growing objects with radii in the 
range RtoR + dR is equal to (j)(L) dL, where L and dL are the luminosity and range corresponding to this interval in R. 
Thus, 

■Wgrow (R) = N b 0(L) f^J . (5) 

Similarly, the distribution of stalled shells is given by: 
N BtM (R) = N b </>(L) (^f\ , (6) 




2.4 Continuous Creation, Luminosity Spectrum 

2.4-1 Growing Shells 

The form of N(R) in this case can be obtained from the generalisation of either equation ^ or ^ In the former case, N(R) 
is obtained by integrating equation ^ over the luminosity distribution: 

AW^) = J i> 4>{L) fej dL , (7) 
whilst in the latter case N(R) is obtained by integrating equation ^ over a continuous creation rate: 

N SIaw (R) = £ ^^ L )(^J dt ■ (8) 

Equations ^ and ^| are equivalent in this situation, provided the limits of integration are correctly chosen. Note that these 
limits are functions of R, determined by the limiting conditions required to produce a given R. For a luminosity distribution 
truncated at L m in and L max , the lower limit of integration for equation |^ is 

Li=mm(L{R,t e ),L f {R)*j , U > L min . (9) 

Here, L(R, t e ) represents the luminosity for which the shell attains the given radius R within the maximum allowable time, 
t e , whereas Lf(R) is the luminosity such that a shell of luminosity Lf stalls at R, i.e., Rf(Lf) — R. Likewise, the upper limit 
of integration is 

L u = L max , (10) 

representing the luminosity of the youngest shells that have grown to radius R. In the same vein, the lower limit of integration 
for equation ^ is 

ti=t{R,L m&K ) , (11) 

implying that no shells contribute to N(R) if they were formed so recently that even the most luminous have not yet grown 
to R. Likewise, the upper limit is 

t u = mm(tf(L), t(R,Lmin)j , t < t e . (12) 
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For the case of continuous growth with no possibility of stalling, the most distant epoch contributing to N(R) is determined 
only by L m i n . Thus in that case t u = t(R,L m i n ), as long as this is shorter than t e . If, however, the possibility of stalling is 
included, then the contribution from growing shells is additionally limited to periods shorter than tf(R,L), where tf is the 
age at which a shell stalls at radius R (see equation ^). 

2.4.2 Stalled Shells 

By analogy to equation |^, the stalled shells are described by, 

N stM (R)= /%<H£) (^A dt , (13) 



dL 

with t t = tf(R), and t u = t e , where tf < t 



2.4-3 Surviving Shells 

In this section, we also include the possibility of a survival phase (§2.1) beyond t e , of duration t s . For shells stalling at tf < t e 



N m (R) = 



i> dt 4>(L) 



dR f 
~dL 



while for those that continue growing until t e , therefore having tf > t a 

/ \ -i 



N sur (R) 
These yield, 



te+*s 



V dt 4>{L) || 



N snr (R) = %bt a 4>{L) I ^ 



and 



N SUI (R) =il)t 3 (j>{L) 



OR 
dL 



(14) 



(15) 



(16) 



(17) 



2.4-4 Total Size Distributions 

The complete shell evolution is summarized by 

R{L,t) , t < t f (L) , t f {L) < te 



R(L,t) = < 



Rf(L) 
R(L,t e 



t f (L) <t <t e +t s 



tf(L)<t e 



t e <t <t e +t B , tf(L) > t e 



(18) 



As mentioned above, any given L uniquely determines corresponding stall parameters Rf(L) and tf(L). Applying this 
relation to t e , we denote the stall parameters corresponding to t e as R e and L e . In other words, a superbubble with input 
power L e would just stall at radius R e and age t e . Since the input power ceases at this time, t e is the maximum possible 
stalling age. Thus shells larger than R e are all either growing or in the survival stage, because their tf > t e . Conversely, the 
population of shells smaller than R e contains a combination of growing, stalled, and surviving objects. 

Therefore, for R < R e , the total distribution is given by, 

N(R) = N gIOV {R) + N stall {R) + N SUT (R) . (19) 
For the case L m m < L(R,t e ) and L ma x — ► oo, the contribution to N(R) from growing superbubbles is given by equation [| 



-ZVgrow (R) — 



t f (R) 



dt 



(20) 



whilst the contribution from stalled systems is given by equation [13f 
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iV s taii (R) =i>-(u- t f (R)) </>(L) fet^ . (21) 

The distribution for surviving shells, N sur (R), is given by equation |l6|. 
For R > R e , the total size distribution is, 

N(R) =N gIOV {R)+N SUI {R) . (22) 

For the contribution due to growing superbubbles, equation |§| is now written: 

iV gr ow(i?) = \ <f>(L) dt , (23) 

again for Z/ m in < L(R, t e ) and L max — > oo. Likewise, the shells in the survival phase are, by assumption, preserved at their 
instantaneous radii at time t e , so that N SUI (R) is given by equation 



3 EXPRESSIONS SPECIFIC TO THE STANDARD EVOLUTION 

We now incorporate the standard evolution for adiabatic, pressure-driven superbubbles to the expressions developed above. 
The growth of the shell radius is described by (e.g., Weaver etal. 1977): 



250 



1/5 



£l """" 5,s/! • (24 » 

where p is the mass density of a uniform ambient medium, and all units are cgs. The corresponding evolution of the interior 
pressure is 

Pi= 3/5 -4/5 (25) 

(3850tt) 2 / 5 1 v ; 

In what follows, it is convenient to use R e , t e , and L e as scaling parameters, recalling that R e and L e correspond to the stall 
criterion at the characteristic time t e . Equations E4 and E5] can now be written, 



R ( L \ / \ 
~RZ = IT \U ) 



/ \ 2/5 , s -4/5 
Pi L \ t \ 



(26) 



(27) 



As discussed previously, we consider that the shell growth stalls when Pi — Po, the ambient pressure. This condition 
therefore implies a stall age: 

1/2 

(28) 



te \ Le 



Thus, by equation |26|, 

Re \ L e I 



(29) 



Equations ^8] and ^ demonstrate the correspondence between stalling age i/ and radius Rf for a given input power L. The 
corresponding relation between tf and Rf is, 



t L = R L 

te Re 

Thus the final age tf is directly proportional to the final radius Rf. 
For a given t e , equations ^i] and [25] give: 

-1/2 



(30) 



Re = (^m H ) n~ 1/2 P 1/2 t e (31) 

550-7T \-3/2 -3/2 n 5/2 ,2 

L e =^-(Afm H ) 'n P ' t e , (32) 
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in cgs units, where n is the number density of the uniform ambient medium and /i is the mean particle weight. Note that 
equations |3l] and are applicable in general for the relation between Rj, tf, and Lf as well. For ISM parameters n = 
0.5 cm -3 , (i = 1.25, and Po — 3 x 10~ 12 dyne cm -2 , the characteristic time t e = 40 Myr implies R e = 1300 pc and 
L e = 2.2 x 10 39 ergs _1 . 

3.1 Continuous Creation, Single Luminosity 

As described in §2.2, the growing shells are the interesting case here, for which equations and ^ imply, 
AT grow (7?) = | ^ |H ^ | • (33) 

Thus for a single- valued LF of luminosity Lo, the size distribution is an increasing function of R 2 ^ 3 , which results from the 
slower growth of large shells. As the stalled shells accumulate, they also superimpose a 5-function at Rf(Lo), described by 
equation ^. 

3.2 Single Burst, Luminosity Spectrum 

Here we assume that 4>(L) is a power law as given by equation j|. For the growing shells with age t < t e , equations ^| and p6| 
yield, 

i-/3 / \ 4-5/3 / \ ~ 3+3 ' 3 

N gtaw (R) = 5AN b (l-F at ) £ j (l) , (34) 

where 1 — F at is the fraction of N b corresponding to growing shells. The fraction of stalled shells is given by, 

F st = J AL~' 3 dL , (35) 

where L is the luminosity corresponding to the largest stalled shells. As we shall see in §5, realistic values of the LF index fall 
in the range 1 & P ^ 3, for which N(R) oc R 4 ~ 5f3 would be an inverse power law in R. This results from the larger fraction of 
small, weak-L shells. 

The stalled shells are described by equation ||. Together with equation this gives, 
-a ( \ 1-2/3 

iV 8 taii(-R) = 2AN b F st j j . (36) 

The total size distribution at times before all the objects have stalled and R < R e would be the sum of equations and 
|3rj . Owing to the large numbers of weak-L objects, equation ^ will normally dominate N(R) in relative numbers. However, 
these stalled shells will only be present out to radius Rf(tb) corresponding to the age if, of the burst. At this radius, there 
will be a discontinuous jump by a factor of |, resulting from the different coefficients of ^ (cf. equation^) and (cf. 
equation £91) . The subset of objects with R > R e will not stall in t < t e and will therefore follow the relation given by 
equation B4l 

3.3 Continuous Creation, Luminosity Spectrum 

3.3.1 R< R e 

We shall now evaluate the expressions in §2.4.4, for the population of shells with radii R < R e , incorporating 4>{L) as given by 
equation |l]. For > |, L m i n < L(R,t e ), and L ma x — * oo, the size distribution for the growing shells is given by equation |2o[ 
With the aid of equation |2(] this yields, 

rl-B . f t>\ a -W f f \ - 2 +3/3 

«=^t^ (I) \i) • ( 37 ) 

For each R we are summing over all ages up to tf that yield a stall radius R, so that (from equation ^) equation |3t] becomes, 

T 1-/3 f / D \ 2 - 2f3 

^ fl )=^-sr^UJ • (38) 

In the case that < |, equation |^ is dominated by ti (equation |ll| ), and if Lmin > L(R,t e ), then t u — t(R, L m i n ). For 
both of these cases, 
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row 

(7?) = 5Atp 



rJ-B 



-t e 



R 



-2 + 3/3 \R 



2/3 



(39) 



Thus for I/ m in 7^ 0, the slope of the size distribution turns over from 4 — 5/3 in equation ^ to a positive slope of |, at 
R < R(L-max,t e ), where R(L m in,t e ) corresponds to the case in the absence of stalling. Beware that the resulting peak in 

N t 

shall see below 



grow(-R) does not correspond to the maximum peak in the total N(R), since the latter is dominated by stalled shells, as we 



For the stalled shells, equations ^ and |2{] yield, 

1-2/3 

' 1 



7-1-/3 

N stM (R) = 2Ai> -4r- t 



Re 



R_ 

Re 



(40) 



Note that for these stalled objects, the term ^1 — jS-^ = ^1 — ~t~^j by equation |3(j. This therefore represents the fraction 
of their lifetime t e that is spent in the stalled state, for shells of radius R. Although multiplying out equation ^ yields two 
terms with dependences _R 1_2/3 and R 2 ~ 2/3 , N sta ,n(R) is dominated by i? 1_2/3 since we are in the regime -jp- < 1. This is the 
same dependence as the stalled objects in equation 

The shells in the survival phase are given by equation which becomes 



7-1-/3 

N SUI {R) = 2Aip -§— t. 

Re 



1-2/3 



(41) 



Again this is similar to equation since this case is equivalent to the distribution resulting from a single burst containing 
ipt s objects. 

Adding together equations |3f| and ^l], the total size distribution for superbubbles with R < R e and /3 > | is: 



N(R) = Aip 



L l-f3 
Re 



1-2/3 



n i \ 9-6/3 



(42) 



3.3.2 R > Re 



For shells larger than the characteristic radius R e , the distribution of growing objects is given by equation G3l This is identical 



to equation |2C], except that the integration is over the range (0, t e 
yields, 

-5/3 



A r grow(fl) = 5Atp 



rl-/3 

"rT 



te 



-2 + 3/3 



rather than (0, t/(L)). For (3 > |, equation |23| therefore 



(43) 



Note that the R dependence of N gIOW (R) in equation [43] is different from that of shells that grow with the possibility of 
stalling (equation ^) . This is due to the fact that the upper limit of integration in equation |^ is a function of R only when 
the possibility of stalling is included. 

The shells in the survival stage are described by equation |l^, yielding, 

-5/3 

AUr(-R) = t s ' 



Re 



Adding together equations ^ and |S[ the overall size distribution for supergiant shells with R > Re and (3 



N(R) = 5Ai/j 



Ll: 



Re 



4-5/3 



-2 + 3/3 



(44) 



(45) 



This again exhibits the R 4 ~ 5/3 dependence as seen in equation which describes the distribution for freely growing shells. 
This reflects the fact that objects in this size range will never attain the stalling criterion. In general, this relation is a 
very steep function of R compared to equation ^ for R < R e . The criterion R > R e corresponds to L > L e , which may 
provide a convenient defining criterion for these larger-scale phenomena, for example, starburst events. There might possibly 
be applications for equation ^5| for a large sample of starburst phenomena with uniform properties. 

In the case that (3 < |, equation |§] is dominated by the lower limit ti, as in the regime for R < R e . N gvov ,(R) in that 
case is therefore given by equation^. Note that, for t s > 0, N(R) is not continuous at R — R e . This is a consequence of the 
change in slope of the L — R relation for shells surviving at Re, as is apparent from equations |26| and §§ (see §2.4.3). 
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3.4 Robustness of the Result 

The dependence N(R) oc R}~ 2 ^ appears to be a fairly robust description, since stalled superbubbles quickly dominate most 
situations for both constant and instantaneous superbubble formation. In addition, we also mentioned above that the slope 
will not change if the breakup of shells into subunits is independent of the original size (Clarke 1996). 

Furthermore, the P 1_2/3 dependence turns out not to be unique to the standard evolution. Consider the general relation, 

(46) 

Within factors of unity, the stalling criterion is equivalent to setting the shell ram pressure p s v 2 ~ Po, where p s is the shell 
density and v is the shell expansion velocity. Therefore, this is the same as requiring that v stall at a limiting velocity Vo 
associated with Po of the ambient medium. We have, 

Setting this equal to vo, we obtain, 

which is the generalization of equation Note that R e is defined such that 22L . = 1. With equation we then have, 
^L = JL.^L (*L\ (49) 

te VO t e \R e J ' V ' 

analagous to equation |3(], The relation tf oc Rf therefore applies generally to any growth law described by equation In the 
standard adiabatic model, x = i and y = |, yielding tf oc L 1 ^ 2 . Other values of x and y yielding = | will also produce 
equations ^ - |3o[ and therefore our previous expressions for iV s t a ii(P) (equation |4o|) and A r sur (P) (equation ^) would still 
be valid. 



3.4-1 Momentum-Conserving Evolution 



It turns out that the non-adiabatic, momentum-conserving shell evolution described by Steigman, Strittmatter, & Williams 
(1975) fulfills this condition. The momentum- conserving growth is given by, 



R_ 

Re 



1/4 , x 1/2 

i) 



(50) 



Therefore, the only component of N(R) that differs from the adiabatic case is that for A r grow (P), for which equations ^ and 
|(] give: 

2-2/3 



iVgrow(P) =4AV> 



Re 



R 



-1 + 2(3 \ R e 



(51) 



in contrast to equation Q The total size distribution for the momentum-conserving evolution in the regime R < R e , (3 > | 
is therefore, 

1-2/3 



N(R) = Aif) 



L, 



1-/3 



/ i' 



2(t e+ts ) + ^Lte 



(52) 



which may be compared to equation |42| N(R) in this case (equation ^) is still dominated by the term having P 1_2/3 
dependence. 

To check whether typical observed superbubbles are still in the regime R < R e for this case, we need to compute the new 
value of R e , again for t e = 40 Myr. The unsealed relations for the radius and interior pressure are: 



R = 



Pi = 



3L 



1/4 



,1/2 



4ttP 2 ' 
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where M and Voo are the mass-loss rate and velocity of the "wind," thereby giving L = 4Mv^,. The pressure equilibrium 
stalling criterion now yields, 

R e = 6 1/2 (umn) ^ n~ 1/2 P 1/2 t e ■ (55) 

The dependences are the same as equation [n], but with a coefficent of 6 1//2 in the momentum-conserving case versus 5/7 1//2 in 
the adiabatic case. Thus, for a given stall age tf, the corresponding stall radius of a momentum-conserving superbubble will 
be larger than that of an adiabatic superbubble. This results from the fact that the same value of tf corresponds to a larger 
value of L in the momentum-conserving case because more energy is needed to maintain growth during the entire period tf. 
Therefore the final stall radius Rf for the momentum-conserving evolution is larger than the adiabatic case for a given tf. 
For the same default ISM parameters as the adiabatic case, equation ^ yields R e = 1700 pc. Most observed superbubbles in 
galaxies indeed have radii smaller than this value. The form N(R) oc 7? 1-2 ^ is therefore a fairly robust result for this simple 
representation of superbubble evolution. 



4 DERIVING THE MECHANICAL LF FROM THE H II REGION LF 

We must now consider the form of <j>(L), the mechanical LF (MLF; equation [l]) , with regard to its slope j3 and the upper and 
lower limits. The MLF results from the stellar census present in the OB associations, which in turn manifests its presence 
through the HII region LF (HII LF). As recognized by Heiles (1990), the MLF and HII LF are closely linked. We will focus 
on clusters that are sufficiently rich so that statistical fluctuations in stellar membership have a negligible effect on the mean 
Ha luminosity per star. We term such clusters "saturated," and consider this limit for now, deferring a discussion of smaller, 
"unsaturated" clusters until §4.3. 

For the saturated clusters, the total initial Ha luminosity Co, and the mechanical luminosity L, are both proportional to 
the total number of stars in the cluster, thus implying L oc Co. Therefore, the slope of the MLF is equal to that of the HII 
LF. This is not, however, necessarily the same as the slope of the observed HII LF (von Hippel and Bothun 1990), because the 
HII region luminosity fades with time as the ionizing stars expire, and the observed HII LF contains objects of differing ages. 
Therefore, before we simply take the empirical slope of the H II LF as that of the MLF, we must first consider the possible 
effect of evolution on the HII LF. We will now examine this effect analytically. 



4.1 Analytic Evolution of the HII Region LF 

We characterise the evolution of the HII region luminosity by 

C = C f(t) (56) 

where /(0) = 1. We thus assume that the fractional fading in a given time, and hence the IMF, is the same for all clusters. 
Since the observed HII LFs are well described as power laws (e.g., Kennicutt, Edgar, & Hodge 1989), we will presume the 
initial HII LF to have the same form: 

<^ = H C 0) =A C -e . (57) 

We adopt the index f3 as in the previous sections, because as argued above, we assume that 4>(L) has the same exponent as 
$(£o). As with (f>(L) (equation [l]), <&(£) is defined such that $(£) dC is the fraction of objects with Ha luminosities in the 
range C to C + dC. If the magnitude of M? is large at early times, we would expect the steady-state observed HII LF to be 
steeper than the initial, since the brightest objects are quickly diminished. On the other hand, if / evolves strongly at late 
times, then the initial slope /3 is largely preserved. 

The evolution of <&(£) may be considered analogous to that of N(R) (equation above: 

$(£)= / £0 UP ^ $(£o) (-^] dC , (58) 



c 



where ip is the constant cluster formation rate as before and the partial derivative with respect to t is evaluated at constant 
C. The observed HII LF is determined by the contribution of objects at all ages to a given luminosity bin, and these therefore 
fall in the range (C, £o,up). From equation |H(| 

■dt =C °di ■ (59) 
So in general, 

<&(£) = / °' UP V *(£//) (-Co dCo (60) 
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where / and Mr are functions of C and Co as described by equation |56[ The upper limit of integration is determined by the 



upper limit to the initial HIT LF, £o,up- For a power law LF given by equation 



$(£) = -A £~ 



df 



(61) 



We now consider the form of f(t). This depends on the stellar mass (m) vs. main-sequence lifetime (t ms ) relation, as 
well as the relative contributions to total C from stars of different masses. Leitherer (1990, Figure 8) shows that fortuitously, 
for a standard Salpeter (1955) IMF, the contribution to C from each unit mass bin 4^-, is approximately constant for all 
m > 20 Mq , while for lower-mass stars it plummets steeply. 

As an exploratory case, we therefore make the approximation that only stars with m > 20 Mq contribute to C, and that 
the different masses contribute equally. We represent the m — i ms relation for m > 20 Mq as a power law: 

. -d 



(62) 



The Geneva stellar evolutionary models (Schaerer etal. 1993; Maeder 1990; Maeder & Meynet 1988) show that this relation 
actually increases somewhat faster than a simple power law, but for our purposes, the approximation is adequate. Fits to 
their models, which are the same as those used in Leitherer's (1990) study, yield d ~ 0.65 in the range 20 < m < 120 Mq. If 
we furthermore assume that stars in each mass bin contribute a constant luminosity for the duration of their main sequence 
lifetimes and nothing thereafter, we then obtain the fading function, 



(t/t 2 Q)- 1/d -l 
(ti/t20)- 1/d -l ' 



t < U 



t > ti 



(63) 



where £20 = t ms (20 Mq), and ti is the time at which the L fading begins, which is i ms (m up ), where m up is the upper-mass 
limit of the IMF. Applying equation |r33| to equation ^l| we obtain: 



$(£) = A ip £~ 



^ -1 f + 1 



- df +U 



(64) 



The last term in equation pj corresponds to shells with 



using a to denote the slope of the observed HII LF, $(£) oc CT 
ages t < ti, for which equation ^ diverges. 

Equation ^ demonstrates that for /3 > d + 1, the integral is always dominated by its upper limit, so that ageing does 
not affect the form of the HII LF in this case. For /3 < d + 1, if is sufficiently large, then $(£) will asymptotically 
approach a slope of d + 1. Values for the observed slope are found to be a ~ 2 ± 0.5 (e.g., Kennicutt etal. 1989). Figure [j] 
shows the resultant <&(£) from equation numerically integrated for initial slopes (3 ranging from 0.5 to 2.5, with d = 0.65, 
-Co, up = 10 40 ergs -1 , and m up = 120 Mq. The curves are normalized to the same upper value to facilitate comparison. Each 
curve is labeled with the input value of (3, which is within 4% of the resultant slope a, fitted over the region log L < 36.5. Thus 
it is apparent that evolution has essentially no effect on these slopes for m up = 120 Mq. With the dashed line, we also show 
the result for /3 — 1, d — 0.65 and m up = 10 7 Mq; this model confirms that the observed slope does tend to the asymptotic 
limit of d + 1 = 1.65 if m up is sufficiently large, although this effect is clearly negligible for realistic situations. 

Alternatively, we can estimate f(t) directly from the population synthesis work of Leitherer & Heckman (1995). Figure 37 
of that work shows the evolution of the total ionizing photon emission rate ^(H ), and hence L, for clusters with different IMFs 
and m up . The fading function is relatively insensitive to these variables and is well approximated by a power law decrease 
with time, which does not drop to zero at £20, but persists with the same power law at late times. From the Leitherer & 
Heckman (1995) standard model with a Salpeter IMF of slope 2.35 and m up = 100 Mq, we fit a fading function slope 77 = 5.0 
in the regime log(t/Myr) > 6.5 and log^H^/s- 1 ) > 44. This slope was also found by Beltrametti etal. (1982). 

We thus parameterize the evolution as: 



t < U 



f 



{t/U)~ v , t>U 
so that integration of equation | 



<£(£) = 



1) 



yields: 

C 

£o,ui 



+ i>uA r 13 



(65) 



(66) 



This again implies that the slope of <!>(£) remains unaffected by ageing for /3 > - + 1, whereas if /? < - + 1, ageing produces 
an observed slope a = 1 + independent of (3 and m up . This minimum value of a reflects our argument above, that the 
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Figure 1. HII region LF resulting from the fading function given by equation |6J, for log(£ up / erg s 1 ) = 40, d = 0.65, and m up = 
120 Mq (solid lines) and m up = 10 7 Mq (dashed line). The input values of /3 are shown for each model. 

observed slope of $(£) is steepened if the nebular luminosity evolution is strong enough. Figure |^ demonstrates this behavior 
for r] = 5 and a range of (3, with the curves again normalized to the same value of £o,up- Thus if the slope of $(£o) is steep, 
intrinsically luminous systems are rare and therefore the faded remnants of such objects make a relatively small contribution 
to any bin of observed £. Note that this steepening behavior in a was not apparent, for plausible values of m up , in the case 
where f(t) was set to at time t = t2o (equation ^). In that case, the dynamic range of the fading function f(t) is small over 
most of the lifetime of an H II region. The stronger evolutionary effect resulting from equation |66| stems from the prolonged 
fading at late times in this prescription. 

Note that 77 essentially corresponds to 1/d in the previous case. Recalling that d is the power law slope for the m — t ms 
relation (equation |u^ ), a value of 77 = 5.0 suggests d = 0.2. This is considerably shallower than our previous estimate of 0.65, 
which was fitted from the model results of the Geneva group {e.g., Maeder 1990), the same stellar models used by Leitherer & 
Heckman (1995). This would suggest that the luminosity evolution of the HII regions is affected by additional factors besides 
the turnoff of stars from the main-sequence, for example, spectral evolution of the stars while they are still in the H-burning 
phase. Clearer understanding of the behavior of f(t) is needed to confirm the applicability of our analysis. 



4.2 Discussion of Evolutionary Effects 

We have shown how in principle ageing can cause the measured H II LF to be steeper than the initial, as cluster fading shifts 
an excess of objects into low luminosity bins. This effect is only manifest for flatter initial $(£0) with low j3, where any bin of 
objects at given current luminosity C contains an important contribution from initially luminous objects formed in the distant 
past. Thus for that case, it is the fading function f(t) that determines the observed HII LF. Conversely, if the initial $(£0) 
is steep, then any bin of current luminosity is dominated by recently formed objects. In this case, ageing is unimportant, and 
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Figure 2. HII region LF resulting from the fading function given by equation |6Jj, for log(£/ erg s 1 )=40, 7j = 5.0, and m up = 120 Mq. 
The input values of f3 are shown for each model, along with the resulting slope a, fitted to the region log£ < 36.5. 



the observed "£(£) has the same slope as the initial. Quantitatively, if the nebula fades according to a simple power law t 17 , 
then 



a — max 



= 1+ - 
V 



(67) 



Thus for example, if T] = 1, then measured values of a — 2 would imply only that /3 < 2. Hence, for any simple power-law 
f(t) that does not steepen substantially at late times, there is a minimum value a m i n , independent of /3, given by equation |u7j 
Therefore, if different galaxies have populations of H II regions characterised by different values of f3, then this variety would 
be manifest only down to a limiting value set by the f(t). An observed value of a mln = - + 1 can therefore be used to 
determine the actual slope r\ of the fading function. For example, a m i n = 1.4 would imply an adjustment in the value of r\ 
from 5 to 2.5, which would require a much more gradual decline in Ha luminosity for the HII regions than is predicted by 
the stellar models. However, r\ is extremely sensitive to the adopted value of a m j n , so this test is unlikely to be practical in 
the near future. Furthermore, our ability to observe a m ln depends on the actual existence of galaxies with HII LFs having 
initial /3 < a m j n . At present the available data on HII LFs in different galaxies suggests a <; 1.3 — 1.4 (e.g., Banfi etal. 1993), 
but there are not yet enough data to evaluate whether there is a lower cutoff in the values of a. This is evidently an effect 
that can be sought in the future as H II LFs become available for a larger sample of galaxies. 

We find, however, that this evolutionary effect of a minimum a is unlikely to be significant for observed HII LFs with 
spectral indices of a ~ 2, if we use f(t) based on plausible population synthesis arguments. We have experimented with two 
prescriptions: one assuming an instantaneous cutoff in £ at t > £20; and the other using continuous fading, f(t) oc t~ v , based 
on population synthesis models by Leitherer and Heckman (1995). In each case, the power-law fading is too steep for the 
measured values of a to be explicable as an ageing effect, and thus we would conclude that a is a direct measure of j3. We 
caution, however, that this result does depend on the assumed IMF and modeled stellar parameters that are used in the 
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Table 1. Observed Parameters 



Galaxy 


a 




HII Ref. a 




h(pc) 


HI Ref. a 


SMC 


1.9 


93 


(1) 


501 


(2000) b 


(3) 


Holm II 


1.4 


67 


(2) 


51 


625 


(4) 


M31 


2.1 


207 


(1) 


110 


120 


(5) 


M33 


2.0 


257 


(1) 


118 


(100) b 


(6) 



a References: 

(1) Kennicutt, Edgar, & Hodge (1989) 

(2) Hodge, Strobel, & Kennicutt (1994) 

(3) Staveley-Smith etal. (1996) 

(4) Puchc etal. (1992) 

(5) Brinks & Bajaja (1986) 

(6) Dcul & den Hartog (1990) 

b Value for the SMC is based on Staveley-Smith's etal. (1996) 
estimate for the SMC HI morphology; value for M33 is as- 
sumed, not measured. 

population synthesis models. Some of these parameters, for example, stellar ionizing fluxes (especially those contributed by B 
stars), are rather poorly known. Furthermore, this determination of /3 from a is based on stars with m <; 20 Mq, which is the 
population ionizing the HII regions. We are thereby assuming that the slope /3 extends to lower-mass stars, which dominate 
the MLF through their SNe. This is based on assumptions about the form of the IMF, and also the constancy of SN power 
for stars of different mass. 

4.3 Unsaturated clusters 

In the limit of small cluster membership number N*, there is a wide dispersion in average Ha contribution per star, according 
to the precise stellar membership of the cluster. As a result, a power law slope in the IMF produces a $(£o) of the same initial 
slope at high Co, but which flattens at low Co (McKee & Williams 1996). This is produced by the behavior of the scatter, 
due to small number statistics, in contributions to Co from a bin of given N*. These contributions scatter symmetrically on a 
linear scale of Co, but logarithmically, they have a larger displacement to smaller values of log(£o), producing a consequent 
flattening of the luminosity spectrum. Monte Carlo simulations by McKee & Williams (1996) confirm that this flattening 
occurs at Co ~ £ mup , the Ha luminosity due to a single star of m — m up . Clusters with Co < £m up are not rich enough to 
saturate the IMF up to this maximum, and thus have varying mean Ha contribution per star; whereas those with Co > £m up , 
again for a constant IMF, will have C tx N*, and thereby a constant Ha contribution per star. Note that we do not expect 
a similar turnover in the MLF because the relative contribution to L of each star is dominated by the SN power, which we 
assume to be independent of m. The correspondence between a and /3 therefore applies in the regime where both Co and L 
scale with the number of stars in the associations. 

For some galaxies, the slope of the HII LF shows a turnover to a shallower value at lower C (e.g., Kennicutt etal. 1989), 
which could possibly be due to this effect, with objects falling below cluster saturation. The value of £m up , corresponding to 
the saturation turnover, is currently estimated at log£ mup ~ 38.0 — 38.5 (Vacca etal. 1996; Panagia 1973). If the observed 
flattening of the HII LFs is caused by this effect, it is important to use the slope a derived from the upper end of the HII 
LF in estimating /3. However, although the observed value of a might cause an underestimate of f3 when fitted over the entire 
range of C, we will see below that, where the measured N(R) is discrepant from that predicted, it requires shallower slopes 
in the MLF than those inferred from the HII LF. 



5 COMPARISON WITH OBSERVATIONS 

We can now compare the prediction with observations for individual galaxies that have been mapped in HI at sufficient 
resolution, and that have measured HII LFs. Table [j] lists observed parameters for nearby galaxies that we have examined. 
Values for the observed slope a of the H II LF, fitted from data in the listed reference, are given in Table |l], along with the total 
number of detected HII regions N Hn having C > 1 x 10 37 ergs _1 , and total number of detected HI holes iV H i- The limiting 
C is that at which the surveys are complete for all the galaxies. We assume that all the H I holes correspond to superbubbles. 
Table hi also lists the assumed HI scale height h, and reference for the HI mapping. 
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Table 2. Inferred Parameters 



Galaxy 


iVtot 


^(Myr" 1 ) 


log A 


/3 a 


/3o b 


<t(/3o) 




a a 


SMC 


1.9 x 10 3 


43 


32.26 


1.9 


1.9 


0.3 


2.8 


2.7 


Holm II 


1.4 x 10 3 


31 


13.98 


1.4 


1.6 


1.2 


1.8 


2.1 


M31 


4.3 x 10 3 


96 


39.53 


2.1 


1.8 


0.7 


3.2 


2.6 


M33 


5.4 x 10 3 


120 


35.90 


2.0 


1.6 


0.1 


3.0 


2.2 



Uncertainties f) ± 0.2. 
b /3o = (<* + l)/2. 
c a p = -1 + 2/3. 



5.1 Slope of N(R) 

Observed slopes for the HII LFs typically fall in the range 1.5 a ^ 2.5 (e.g., Kennicutt etal. 1989), and thus imply a similar 
range in /3. We therefore see that we are indeed in the regime /3 > | for which we developed the solutions of N(R) above. 
Focusing again on the superbubble size distribution for R < R e (equation we find that such values of /3 imply an N(R) 
that is dominated by the term with a dependence of J? 1-2 ' 3 . This term describes stalled shells, including those that are in the 
survival phase after stalling. Our assumed range of /3 therefore implies a range in slope a, where N(R) oc R~ a , of 2 ^ a Zz, 4. 

Figures ^| and |^ show the histograms of HI hole radii in pc, with Figure ^| binned and displayed linearly, and Figure ^ 
binned and displayed logarithmically. The data in the figures are superimposed with a dashed line showing a least-squares, 
power-law fit to the data, of slope a . This fit was derived from the logarithmic binning in Figure [l|, weighted by the inverse 
of the root-iV errors shown in the figure. Bins without error bars are deemed incomplete and were not included in the fit. 
We also show the prediction of equation [li] in the solid line, having a slope of essentially a p = 1 — 2/3. Where known, the 
HI scale height h is indicated with the vertical dotted line, and the HI survey resolution limit is indicated by the vertical 
dashed line. Values for ce p are given in Table [| and for convenient comparison, a is listed there as well. We also give values 
for /3 = (a + l)/2, predicted from the observed a , and the uncertainties cr(/3 ) derived from the formal standard deviation 
on the fit to a . 

The comparison between predicted and observed slopes is surprisingly good, and is in agreement within the errors for all 
the galaxies. The SMC has by far the most reliable and complete HI data, and shows a v = 2.8 ± 0.4, in excellent agreement 
with ot — 2.7 ± 0.6. This is more vividly demonstrated by the comparison that /3 = 1.9 ± 0.2, taken from the HII LF, 
which is in superb agreement with /3 = 1.9 ± 0.3, predicted by the HI hole size distribution. The agreement for the SMC is 
extremely encouraging, given the high level of completeness in the data for this galaxy (see below). While the agreements for 
Holmberg II, M31, and M33 are additionally gratifying, unfortunately the poor statistics in the data for those galaxies render 
the comparison less significant. 



5.2 Normalization of N(R) 

The normalization of N(R) depends on the coefficient A of the MLF and the superbubble formation rate ip. A can be calculated 
by integrating equation hi over the relevant range of L. This yields, 



Ll-"-L^\ . (68) 



A = (1-/3) 

The lower limit of integration, L m i u , is the power associated with a single SN explosion. In the standard treatment for the 
SN power used above, this implies Z/ m i n = E$\/t e = 8 x 10 35 ergs -1 , where £51 = 10 51 erg, the assumed energy of a single 
SN explosion. In this limit, the approximation breaks down in its treatment of an energy input over the period t e , and we 
will address this point in §5.3 below. However, this does provide a useful value for I/ m i n - We use L e as the upper limit of 
integration because we are comparing the observations to equation^, which applies only for R < R e , or equivalently, L < L e . 
Since N(R) is dominated by stalled objects at small R, and therefore small L, equation |i8] is not sensitive to the choice of 
upper limit. 

The superbubble formation rate %j) = Nt t/(t e + ts), where JVtot is the total number of superbubbles at any given time. 
This assumes ip to be constant over the timescale of t e — 40 Myr, relevant to the observations. As can be seen in Table |l| 
it generally appears that for the galaxies examined, -/V H n > JVhi- This is presumably a selection effect due to the sensitivity 
and resolution limits of the HI surveys. The incompleteness is also better understood for iV H ii than AT HI , so we will therefore 
estimate iVtot from JV H n. Since L for the nebulae diminishes to below the detection limit well before t e , we must apply a 
correction to iVnii to obtain the total number of superbubbles JV to t, which we estimate as follows. 

Unsurprisingly, the proximity of the SMC, surveyed at high sensitivity in HI (Staveley-Smith etal. 1996), reveals a much 
greater surface density of superbubbles than the other galaxies. 7V H i > 2Vhii for the SMC, which, as mentioned above, would 
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Figure 3. Histogram of HI hole radii. Power-law fits to the data with slope a Q are shown with the dashed line, from which inferred 
values of the MLF slope /3 Q are computed. The solid lines show the predicted size distribution with slope ct p , computed with the value 
of fi from the H II LFs. H I scale heights, where known, are indicated by the vertical dotted lines, and the survey resolution limits are 
shown by the vertical long-dashed lines. 
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be expected generally. Furthermore, the excellent agreement in observed and predicted slope for N(R) suggests that for the 
SMC, we can simply normalize the predicted N(R) to the data, obtain N to t, and thereby estimate Ntot/N H n. Assuming that 
the shell census for the SMC is complete for R > 100 pc, we can use the integral / 10 q N(R) dR to normalize the prediction 
for N(R) and obtain JVtot. Integrating equation p2| yields, 



-20 _ t,2-20 
n 100 



(69) 



where iVioo = 159 is the number of SMC shells with R > R 10 o, which is 100 pc. Using /3 = 1.9 from the HII LF for the SMC 
(Tables and ||), this yields N to t = 1.95 x 10 3 , and ip = 43 Myr -1 with t e +t s = 45 Myr. We therefore obtain N to t/N H u = 21.0. 

General application of this conversion factor does assume that N a i is essentially complete at R > Rioo, although Staveley- 
Smith etal. (1996) believe there is still substantial incompleteness in their data. However, based on the following, we expect 
that most of the incompleteness is in the small shells, where the resolution limit of 28 pc becomes important. The SMC ratio 
of N tot /N Hn = 21.0 implies a typical HII region age of 2.1 Myr, for t e + t 3 — 45 Myr. We note that if HII regions fade 
according to the power law prescription (equation ^), the average age of nebulae that are observed before fading below a 
given completeness limit L c is given by, 

_ /°° 4>(L) te(L) dL 

• , (70) 



/~ <P{L) dL 

" (A) 



-i/n 

where t c = U ( -r- ) , so that 



U . (71) 



The mean age i is therefore close to ti because, with the assumed luminosity evolution that fades strongly after ti, most of 
the nebulae will spend only a short fraction of their lifetimes in the fading period before disappearing below the completeness 
limit. For (3 = 2, and taking t ms (m up ) to be ti — 2.8 Myr (e.g., Maeder 1990), equation fn] gives i = 2.2 Myr. Thus the mean 
HII region age based on our assumptions for timescales and luminosity evolution is in excellent agreement with that of 2.1 
Myr implied by the SMC data. We are therefore fairly confident that this ratio of N tot /N Hn may be used as a reasonable 
conversion factor to obtain N tot . However, we caution that applying the SMC conversion universally assumes the same relative 
extinction for nebulae in all the galaxies, which is actually likely to vary. The SMC itself appears to have a large column depth 
of HI (Staveley-Smith etal. 1996), and may have a disproportionate fraction of HII regions lost below the completeness limit 
of the nebular census. 

Using iVtot /Nun = 21.0, we estimate Ntot for the remaining galaxies. These are listed in Table [2], along with the resultant 
tp and MLF parameters A, (3. These parameters can now be used in equation ^ to quantitatively predict the superbubble 
size distributions for each galaxy, that are shown in Figures || and [|. 

The relative normalizations of the observations and predictions are in worse agreement than the slopes, reflecting the 
huge incompleteness in the HI hole data described above. There is also significant uncertainty in the conversion AT to t /Nun, 
which as we argued, may be slightly higher in the SMC than the other galaxies. Hence the disparity in the total numbers of 
predicted and observed N(R) may be slightly exaggerated in these other galaxies. 

One check on the critical parameter ip and resulting normalizations is a calculation of the expected core-collapse (Type 
II + lb) SN rate, ?/;sn = Nsn4>i that is implied by our analysis. The average number of SNe per cluster JVsn, can be estimated 
as 

N m = ±± = J- . (72) 

This again assumes that the stellar winds contribute little to L. The mean logL is set by j3 for constant L m i n . We tabulate 
log L and ^sn for the galaxies in Table |§| with SN rates given in SNu = number of SNe per 10 10 Lq(B) per 100 yr. Table | also 
shows for comparison the empirically estimated core-collapse SN rates ^sN(obs) as reviewed by van den Bergh & Tammann 
(1991). There is a factor of a few uncertainty in these rates, but reassuringly, our computed ipsnifi) are m reasonable agreement 
or less than the empirical estimates. Our normalizations for N(R) are therefore generally quite consistent with the expected 
SN rates, and at least are not substantial overestimates. 



5.3 Peak in N(R) 

The observations show a peak in N(R), which is undoubtedly related to the resolution limit of the HI surveys. However, 
it turns out that the theoretical size distribution peaks at similar radii. Since N(R) is dominated by the stalled shells, the 
minimum stalled R = i? m i n is simply the stall radius corresponding to the minimum input power L m - m , which corresponds 



© 1997 RAS, MNRAS 000, §-?? 



Superbubble Size Distribution 19 



Table 3. Supernova Rates 



Galaxy 


M s (Ref.) a 


UP) 


V>Sn(/3) 


^SN(obs) b 






(ergs" 1 ) 


(SNu) 


(SNu) 


SMC 


-15.7 (1) 


8.6 X 10 36 


1.1 


1.3 - 3.0 


Holm II 


-16.6 (2) 


6.3 x 10 37 


2.5 


(2.9) 


M31 


-21.6 (3) 


4.8 x 10 36 


0.006 


0.03 


M33 


-18.4 (1) 


6.3 x 10 36 


0.19 


0.4 



a References for Mg: 

(1) van den Bergh &: Tammann (1991) 

(2) Puche etal. (1992) 

(3) Sandage & Tammann (1987) 

b From van den Bergh & Tammann (1991); value for Holm- 
berg II is that estimated for galaxies of type Sdm - Im, with 
H = 75 km s" 1 Mpc" 1 . 

to that of individual SNe, as shown above. For R < R m in, N(R) oc R 2 ^ 3 , as given by equation ^9| (§3.3). Our value of 
Lmin = 8x 10 35 ergs -1 thus yields J? m in = 30 pc (equation However, the evolution of individual SNRs is rather different 
from that of superbubbles with continuous energy injection. Our limiting extreme in superbubble evolution yields an SNR 
stall radius that is a factor of 2 or 3 less than estimates by Cioffi & Shull (1991), who consider more realistic SNR evolution 
and slightly different ISM parameters. Consideration of an ambient magnetic field by Slavin & Cox (1992) produces maximum 
radii of ~ 50 pc, 20% less than unmagnetized evolution. So it would seem that the peak in the size distribution of HI holes 
should fall between 30 and 100 pc. 

The HI surveys seem to suggest a peak of ^ 50 pc. The resolution limits have similar values, so at present it is unclear 
whether N(R) continues to increase toward smaller radii, or whether the peak is actually falling in the expected range. It is 
interesting that for the SMC in particular, which has the highest quality data and resolution limit of 28 pc, that N(R) still 
shows no sign of falling off at low values. However, our treatment of superbubble evolution assumes a universal lifetime for 
all shells of t e — 40 Myr. Again, this assumption breaks down for individual SNRs. The peak in the size distribution may be 
affected by their shorter lifetime and shrinkage by the ambient pressure or magnetic tension (Ferriere et al. 1991) after the 
maximum sizes are attained. At any rate, assuming that the peak radius does actually correspond to the typical endstage 
radius for individual SNRs, the observed peak value of R could help constrain SNR evolution, ISM conditions and/or the 
typical SN energy in galaxies for which it is well-determined. 

6 SUPERBUBBLES AND THE ISM 

It is extremely encouraging that the predictions and observations for the slope of N(R) are in agreement for the HI holes in 
these galaxies, especially the largely complete sample for the SMC. Since these superbubbles constitute one of the principal 
forms of structure in the ISM, the implication is that its large-scale structuring is indeed most likely determined largely 
by these OB superbubbles. Furthermore, a strong agreement, as might be the case for the SMC, implies that no additional 
fundamental process is necessary to explain the creation and evolution of observed HI holes. 

The agreement of our relation for N(R) with the SMC data additionally suggests that most of the simplifying assumptions 
made in §2.1, such as constant IMF, single-burst star formation in clusters, uniform ambient medium, etc. are practical in this 
analysis. We have confirmed that the MLF spectrum, in particular, is a critical parameter to N(R), and affects the resulting 
ISM structure. In comparing with Galactic HI observations, Bregman, Kelson, & Ashe (1993) found discrepancies in ISM 
structure functions modeled from a uniform size distribution of HI holes. It would be interesting to see whether accounting 
for a power-law N(R), as required by the MLF, can more accurately reproduce the observations. 

We do caution that Staveley-Smith et al. (1996) note a strong coherence in the dynamical ages of the SMC shells, suggesting 
a short formation burst. This is in contradiction with our assumption of constant ip. We predicted the form of N(R) for a single 
burst in §3.2, finding that N(R) cx T? 1 " 2 * 3 for R < R f (t b ), where t b is the age of the burst. For R > R f (t b ), N(R) oc R 4 ~ 5f3 , 
corresponding to growing shells, with a factor | discontinuous jump. Staveley-Smith etal. (1996) find tb = 5 Myr for a 
burst scenario, yielding Rf — 160 pc. The data in Figure ^ extend to twice this radius, with no suggestion of the predicted 
discontinuity. The observations therefore appear to be inconsistent with this single burst. Since the constant formation model 
implies that most shells have stalled, we suggest that the Staveley-Smith etal. (1996) result may stem from this majority of 
stalled objects having spurious expansion velocities attributed to them. The stalled shell walls may instead reflect random 
ISM velocities. We plan a followup study of the distribution in expansion velocities, which will clarify this issue. 

As is apparent in Figure ^, the two disk galaxies, M31 and M33, show greater disagreement in predicted and observed 
slopes of N(R) than do Holmberg II and the SMC, which are both Magellanic irregulars. Owing to the apparent enormous 
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incompleteness in the HI hole samples and high uncertainties on the fit to N(R), we cannot attribute significance to this pos- 
sible correlation between galaxy morphology and relative agreement in slopes. Nevertheless, the contrast in relative agreement 
between the different galaxy types is suggestive and interesting. Since the gas distribution and dynamical processes in spiral 
galaxies differ from those of dwarf irregulars, it is reasonable to suspect that superbubble evolution is likely to differ between 
these morphological types of galaxies. For example, Ferriere (1995) shows that the three-dimensional evolution of superbubbles 
in an exponential gas disk is quite different from the spatially uniform expansion we have assumed here. Differential rotation 
in the disks and radial effects such as the distribution of gas and star-forming regions are other examples of factors that are 
likely to be important in spiral galaxies. 

6.1 The Smallest and Largest Shells 

The regimes where we expect significant disagreement between the prediction and observations would be at the extremes in 
radius. As described above, the smallest shells should correspond mostly to individual SNRs, whose evolution is significantly 
different from what is assumed by equation The continuous wind approximation is also likely to break down for superbubbles 
created by only a few discrete SNe. The sizes of such objects may be underestimated since observational and theoretical 
evidence suggests that many SNRs strike the superbubble walls, therefore converting their energy into direct kinetic impulses 
and shell X-ray radiation instead of thermal energy in the shell interior (e.g., Oey 1996; Chu & Mac Low 1990; Franco 
etal. 1991). In addition, we may expect a significant contribution from Type la SNe. However, the addition of Type la SNRs 
to N(R) would cause a jump in the peak predicted at i? m i n . There is as yet no evidence for such a jump, again demonstrating 
that J? m in itself is empirically not yet apparent. 

At large R, we expect the HI scale height h to be a critical factor in the evolution of the superbubbles. The growth of 
the superbubble radius to h allows the interior hot gas to break out of the galactic disk and depressurize the shells. Heiles 
(1990) assumed a bimodal evolution, where for R < h, the shell growth followed R <x t 3 '' 5 given by equation [m]; and for 
R > h, the growth was described by a coasting, momentum-conserving phase with R tx t 1 ^ 3 (cylindrical geometry). Such an 
effect should cause a steepening in slope a of the size distribution at large R: since these shells can no longer grow as large 
as they would have done adiabatically, they accumulate at smaller final radii. However, the observations are showing slopes 
that are generally shallower than a p , rather than steeper, thus the limiting effect of h is not apparent in the data for these 
galaxies. As seen in Figure | N(R) does appear rather irregular in M31, especially at R a few tens of pc above h, hinting 
at the possible effect of breakouts. However, for M33 and the SMC, N(R) follows a power-law distribution quite smoothly, 
even for R > h. There is a hint of a break in the SMC data around 300 pc, that could possibly correspond to an equivalent 
h, which is indeterminate for this galaxy. 

It seems likely that some other mechanism might also create superbubbles with R > h, generating an excess of large 
superbubbles, and thereby flattening the observed slope a of the size distribution. As discussed by Heiles (1990), among 
others, the two principal candidate mechanisms are propagating star formation and infalling high- velocity clouds (HVCs). 
Propagating star formation (e.g., McCray & Kafatos 1987) produces generations of OB associations in close spatial proximity, 
thereby increasing the mechanical power L and extending the duration of energy input beyond t e . This can therefore lead to 
supergiant shells with radii larger than are likely to be due to individual OB associations. Impacts by HVCs have also been 
demonstrated to produce shell structures with radii of 10 2 — 10 3 pc (e.g., Tenorio-Tagle etal. 1987; Rand & Stone 1996). 
But these mechanisms, if applicable, apparently do not create enough shells to clearly distinguish their contribution to the 
superbubble size distribution in the current data, and blend smoothly with N(R) predicted by the MLF for OB associations. 

It is also quite possible that many of the larger shells are created by the merging of smaller ones. This would cause a 
flattening in the slope a of the size distribution, since smaller shells would be eliminated to combine larger ones. The inter- 
cluster distance for OB associations is typically a few hundred pc in these galaxies, so we expect merging to be important for 
shells with radii larger than this range. We note that the observed slopes are all flatter than predicted, with the exception 
of that for Holmberg II, which has by far the worst uncertainty. We commented above on the possibility of disagreement 
between predicted and observed slope in the spiral galaxies versus the Magellanic irregulars; another potential explanation 
for such a discrepancy is that the concentration of star forming regions in the spiral arms enhances merging and propagating 
star formation, preferentially encouraging the production of the very large shells and thereby flattening the resultant a. 

6.2 Porosity of the ISM 

The favorable comparison between the observations and prediction encourage us to apply this analysis to the global structure 
of the ISM in galaxies. A fundamental outstanding issue is the relative importance of hot, coronal gas, which presumably 
originates in superbubbles and Type la SNRs, in relation to the cooler phases of the ISM. This question is traditionally 
addressed by means of the porosity parameter Q (Cox & Smith 1974), which is the ratio of total volume or area occupied by 
superbubbles to the total volume or area of the host galaxy. Following directly on the preceding section, we can also evaluate 
the degree of shell merging or overlap by examining Q. 
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Table 4. Porosity Parameters 



Galaxy 


Rg 


Predicted 


Observed 




(kpc) 


Q2D 


Q3D 


Q2D Qw 


SMC 


2 


2.1 


0.3 


1.6 0.1 


Holm II 


7 


1.5 


1.2 


0.1 0.08 


M31 


17 


0.03 


0.06 


0.01 0.02 


M33 


7 


0.3 


0.9 


0.08 0.2 



Following Heiles (1987, 1990), we compute both the two-dimensional Q2D, and volume porosity Qzu- Q20 is especially 
appropriate for the disk distribution of OB associations in spiral galaxies. We have, 



Q2D = A„ 



A b N(A b ) dA b 



(73) 



where At, = irR 2 is the projected superbubble area, and A g is the area of the galactic disk. Rewriting this expression in terms 
of N(R) dR: 

"-Re 



Q2D = {nRg)' 



ttR 2 N{R) dR , 



recalling that we are in the regime R < R e , j3 > | for the standard model. We obtain, 



Q 2 d = Rg 2 AipL 1 - 13 R\ 



B C 
+ ■ 



4-2/3 5-2/3 



Rmin 
Re 



4-2/3 



B C 
+ 



4-2/3 5 — 2/3 V R 



Rmin 



where R g is the radius of the galactic disk, and 

B = 2{t e +t a ) 
9-6/3 



C 



(74) 

(75) 

(76) 
(77) 



-2 + 3/3 

For j3 > 2, Q2D is dominated by Rmin, but for /3 < 2, the largest shells dominate. Interestingly, the relevant values of /3 fall in 
this transition. Thus for some galaxies, the few largest shells dominate Q2D, whereas in others it is dominated by the many 
individual SNRs. For those with ~ 2, the relative superbubble sizes contribute fairly equally in determining Q2D- 

The C term, corresponding to the term in equation [li] with dependence R 2 ~ 2f3 , is small for parameters of interest, and 
will only dominate for /3 < 1. We may therefore approximate, 

-2/31 

(78) 



Q 2D ^ R~ 2 AxpL 1 - R\ 



2(t e +t s 



2/3 



Rn 



Re 



By analogy, we also compute the three-dimensional porosity parameter: 

5-2/3 



3hRj 



AtpLl' 13 Rl 



B 



C 



2/3 



2/3 



Rn 



Re 



B 



C 



2/3 6-2/3 



Rn 



Re 



(79) 



using the HI scale height h as the relevant galactic vertical extent. i? m i n dominates for j3 > 2.5, so Q3D is almost always 
dominated by the largest shells. Again, we find that the C term is usually small, dominating only for /3 1, so we may 
approximate, 

5-2/3" 

(80) 



Q 



3hR 2 g 



AipL, 



1-/3 



R, 



2(te+t s 



5-2/3 



1 - 



Re 



We take R g to be the distance out to which HI holes are detected in these galaxies. The adopted R g , and resultant Q20 
and Q30 are given in Table ^. These are the exact values, computed from equations [r5| and [T9| We also list the observed 
values of Q2D and Qsn, computed from the HI holes with R > fl m j n . We caution that the predicted porosities could be 
overestimates, since equations [75| and [79] assume a size distribution extending to R e = 1300 pc. Because of the small numbers 
of large shells predicted, it is unclear whether the distributions indeed may be considered to extend to R e . 

For M31 and M33, we take Q2D to be the relevant porosity parameter, and for the SMC, Q-jd- The values in Table [| 
of the relevant porosity parameters for these galaxies are all significantly less than 1. We therefore do not expect a great 
deal of superbubble overlap and merging in these galaxies, thereby limiting the extent and networking of the hot gas that 
originates within the shells. On the other hand, the predicted porosity parameters for Holmberg II are 1. This galaxy 
has a very large scale height, so both Q2D and Q30 axe relevant, as is indicated by their similar values. These predict that 
Holmberg II, unlike the other galaxies, does have an ISM that is dominated by coronal gas. This supports the expectation 
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Table 5. Parameters for the Galaxy a 

Model JVtot V log -4 P a p L(JS) -^sn b Q2D Q3D 

(Myr" 1 ) (ergs" 1 ) (SNu) 

MW96 6.5 x 10 3 1.4 x 10 2 35.90 2.0 3.0 6.3 x 10 36 0.05 0.1 0.3 

SK89 3.8 x 10 4 8.5 x 10 2 46.78 2.3 3.6 3.1 x 10 36 0.14 0.2 0.4 

E2.0 4.0 x 10 5 8.9 x 10 3 35.90 2.0 3.0 6.3 x 10 36 (3.0) 6.7 19 

E2.3 8.0 x 10 5 1.8 x 10 4 46.78 2.3 3.6 3.1 x 10 36 (3.0) 3.9 8.1 

"■Assuming R g = 13 kpc, h = 100 pc, and L q ,b = 2.3 x 10 10 . 

b Value for models E2.0 and E2.3 is the empirical estimate of Vsn f° r the Galaxy 
(van den Bergh & Tammann 1991). 



that the importance of the HIM depends on the level of star formation relative to galaxy size and interstellar conditions. It 
is important to keep in mind that, as seen in §5.2, the normalizations of N(R) are uncertain, and these affect the derived 
porosity parameters, as does the assumption of a distribution extending to R e . However, given the general consistency with 
the SN rates, and possible overestimate of the Ntot/Nan conversion factor, we imagine that the porosities in Table ^ should 
not be hugely underestimated. Our derived Q20 nevertheless essentially agree with the values estimated by Heiles (1990) 
for M31 and M33. His analysis is based on somewhat different criteria with greater uncertainties. For example, he examined 
primarily superbubbles that break out of the galactic disk, and included an ad hoc adjustment in his C — N* relation. 

Note that the existence of the MLF spectrum, which accounts for the clustering of core-collapse SNe, has a potentially 
important effect on the porosity. Whereas N t uniformly distributed SNe will, at their final radii, contribute to Qzd simply 
as N*L^? oc 7V», the same SNe concentrated into one cluster will contribute as (iV*_L m in) 3/ ' 2 oc A^ 2 (for N*L m i n < L e \ 
equation ^). Clustered SNe in superbubbles therefore produce a larger Q3D than the same number of individual SNe. 
However, this is not the case for Q2D, where for both individual and clustered SNe, Q2D oc iV*Z/ m m. Although, as we have 
seen above, the largest shells constitute a larger component of Q2D, the relative degree of clustering will not affect the actual 
value of Q2D- Ferriere (1995) also finds a strong effect of clustering on Q30, using a more complex model for shell evolution 
in an exponential gas disk. 

It is also important to bear in mind the contribution of Type la SNRs (e.g., Slavin & Cox 1993), although in most 
star- forming galaxies, these are outnumbered by factors of 3-10 by core-collapse SNe (van den Bergh & Tammann 1991). 
Heiles (1987) emphasizes that the distribution of Type la SNe should vary with galactocentric radius in disk galaxies, owing 
to the distribution of progenitors, and could therefore cause important radial effects in the ISM porosity. This effect could be 
a factor in the discrepancy between observed and predicted slopes for M31 and M33. As mentioned in §6.1, the HI data do 
not yet provide empirical evidence for the contribution of Type la SNRs in the galaxies we have examined. 

It is interesting to apply our analysis to the porosity of the Milky Way. The HII LF for the Galaxy has been compiled 
by McKee & Williams (1996; hereafter MW96) and Smith & Kennicutt (1989; hereafter SK89), where both sets of authors 
used primarily the compiled data of Smith, Biermann, & Mezger (1978). The HII LF estimated by MW96 (their Figure 2), 
shows a slightly natter slope of a = 2.0, compared to a = 2.3 reported by SK89. The data from these papers yield parameters 
for the Galaxy shown in Table ^. The values of (3 from the H II LF predict a fairly steep slope for N(R) in the Galaxy of 3.0 
and 3.6. For the SK89 data, we base our calculations on an estimated total of 145 HII regions having log£ > 37.84 ergs^ 1 
(see their Figure 1). It is apparent that the predicted SN rates V'sn are an order of magnitude smaller than the observed i/>sn 
of ~ 3 SNu estimated for the core-collapse SNe in the Galaxy (van den Bergh & Tammann 1991). We therefore also compute 
implied parameters based on the empirical SN rate, for = 2.0 (model E2.0) and = 2.3 (model E2.3). Resulting porosities 
are computed assuming R g ~ 13 kpc, based on the data of Smith etal. (1978), and h = 100 pc (Kulkarni & Heiles 1987). 

Unfortunately, this discrepancy in the predicted and observed SN rates has a critical impact on the Galactic porosity 
parameters. As seen in Table [B|, the porosity estimates based on the observed Galactic H II LF are consistent with the low 
porosities estimated for the other galaxies (Table |^). However, the Galactic estimates based on the much larger, observed i/>sn 
produce porosities 2> 1, implying a strong dominance of the hot interstellar component. It is unclear how this discrepancy 
should be reconciled. On the one hand, the observed Galactic SN rate is consistent with the presumed Hubble type around Sb 
(van den Bergh & Tammann 1991). On the other hand, van den Bergh finds a similar discrepancy when predicting ipsn from 
the local IMF, which also underestimates the observed ipsN by an order of magnitude (van den Bergh & Tammann 1991). We 
note that Tammann does not find such a discrepancy, and also that MW96 do find consistency between their HII LF and 
the Galactic SN rate. It is interesting to note that the only other Sb galaxy in our study, M31, shows a similar, but smaller, 
discrepancy in observed and predicted tpsN- 

Other authors have estimated the porosity of the Galaxy based on similar analyses. Ferriere (1995) finds Qzu ~ 0.2 based 
on a more complex, 3D model for shell evolution in an exponential gas disk, including an important shell contraction phase. 
Slavin & Cox (1993) estimate a value of Q30 ~ 0.18 for individual SNRs, but exclude the contribution of larger superbubbles, 
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which we found above to dominate the porosity in our analysis. Heiles (1990) estimates Q%t> = 0.30, ignoring the contribution 
of shells that do not break out of the Galactic disk. Since the measurements of j3 for the Galaxy fall in the regime where small 
shells dominate Q20 , this value is also likely to be a substantial underestimate. Our results seem very broadly consistent with 
previous studies, but highlight the uncertainties in Galactic parameters. Apparently we have simply restated the problem of 
whether or not the hot component of the ISM should strongly dominate in the Milky Way, since this argument is based on 
the observed supernova rate in the Galaxy (Cox & Smith 1974; McKee & Ostriker 1977). 



7 CONCLUSION 

We have used the standard, adiabatic shell evolution to predict the differential size distribution N(R) for populations of 
OB superbubbles in a uniform ISM. The results are strongly dependent on the inclusion of a power-law MLF for the OB 
associations. Another fundamental ingredient is the criterion that the shell growth stalls upon pressure equilibrium with 
the environment. This condition, along with the given characteristic time t e , determines the characteristic radius R e , that 
divides the superbubble population into two regimes of solutions corresponding to N(R < R e ) and N(R > R e ). Re and the 
corresponding L e therefore make convenient defining criteria for a different scale phenomenon, for example, a starburst event. 
For the condition of constant shell creation and power-law MLF, in the regime R < R e , N(R) is given by equation which 
has contributions from growing, stalled, and surviving objects. For reasonable values of the MLF slope (3 ~ 2 ± 0.5, the size 
distribution is dominated by stalled objects, yielding N(R) oc _R 1_2/3 . This R 1-213 dependence appears to be fairly robust, 
since it applies to the standard, momentum-conserving evolution (Steigman etal. 1975) as well. It also describes N(R) for 
the single-burst creation of shells, for R smaller than the stall radius associated with the burst age. On the other hand, in 
the regime R > R e , we find the dependence N(R) oc R 4 ~ 5/3 (equation |45[), a much steeper relation essentially composed of 
growing objects. 

To estimate j3, we adopt the observed slope a of the HII LF. However, the observed a could in principle be steeper than 
the initial f3 depending on the relation between the luminosity fading of the H II regions and the initial slope (5. We investigated 
this problem for power-law nebular fading, and found the existence of a minimum slope a m i n = 1 + - (equation ^5?]), which 
is determined by the power-law index r\ of the fading function. Thus, for /3 > a m i n , we may indeed adopt j3 = a, but for 
[3 < a m i n , a provides only an upper limit to /?, as the observed slope is indeed steepened to the value a m - m . Observed slopes of 
galactic HII LFs are all greater than the predicted a m i n , so these should provide a fairly reliable measure for the slope of the 
MLF. It will be interesting to see whether a lower cutoff in the distribution of observed a in galaxies will become manifest at 
Ornin- We also found that the existence of this minimum a m i n depends on a long-term, power-law nebular fading law, whereas 
if the HII region luminosities are instantaneously extinguished at a specified age, no a m - m will be observed. We also derive a 
useful expression for the mean HII region age of a complete population of nebulae brighter than a given luminosity. This is 
given by equation [n], in terms of j3 and rj. 

We compared the predicted size distributions with observed HI hole distributions in four galaxies. The recent data for 
the SMC (Staveley-Smith et al. 1996) appear to be largely complete, and show excellent agreement with our predicted relation 
of N(R) oc J? 1-2 ' 3 . Despite the fact that our critical assumptions about e.g., the endstage evolution, are highly uncertain 
and crudely treated, the slope of the size distribution in this galaxy can be entirely explained by our prediction. No other 
fundamental processes are necessary to explain the observed HI hole distribution in the SMC. This furthermore suggests that 
our assumptions of constant IMF and coeval star formation in OB associations are broadly useful on global scales. 

The observed and predicted slopes of N(R) are also in agreement for the three other galaxies we examined. It is premature 
to draw conclusions based on these other comparisons, since the HI hole data are largely incomplete, as evidenced by the 
relative numbers of H II regions vs. H I holes. The predicted SN rates based on our analysis also support this conclusion, and 
are in reasonable agreement with observed rates. However, it is intriguing that the spiral galaxies appear to suggest greater 
disagreement in the slope of N(R) compared to the Magellanic irregulars. We anticipate that radial and disk properties 
associated with spiral galaxies should cause differences in the observed superbubble size distributions. The effect of the disk 
scale height should be especially apparent in these galaxies, limiting the growth of shells at these sizes, yet such an effect 
is not apparent in the data, which show plenty of shells at larger radii. Processes such as propagating star formation and 
merging could contribute large objects that counteract the expected dearth of large shells. At present, however, the data for 
the spiral galaxies, especially M33, are adequately fit by a simple power-law distribution, and there is no strong evidence for 
a change in characteristics of large superbubbles. 

We also predict a peak in the size distribution corresponding to the smallest stalled shells, which would be individual 
SNRs. The spatial resolution of the data do not yet allow meaningful comparison with prediction. The contribution of Type la 
SNRs, which could be substantial (Heiles 1987), should also be apparent near this peak. 

Our derivation is easily applied to estimating the porosity of the ISM. We find that, not including Type la SNRs, Q-m is 
dominated by individual SNRs for f3 > 2, but the few largest shells for < 2. On the other hand, Q3D is usually dominated 
by the largest superbubbles, rather than the multitude of small shells and SNRs. Our estimates for the porosities in these 
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galaxies generally show values substantially < 1, with the exception of Holmberg II. This therefore predicts that this galaxy 
would be dominated by a hot interstellar medium, whereas the others would not. Furthermore, merging would not be expected 
to be a dominant process in the galaxies with porosities <C 1, although merging could still flatten the slope of N(R) if it is 
encouraged by spatially clustered shell distributions. Our porosity estimates for these external galaxies are in good agreement 
with previous calculations (e.g., Heiles 1990), with the caveat that they are dependent on the normalizations for N(R). 

However, in predicting the SN rate and porosity of the Milky Way, we find a critical discrepancy between the predicted 
and observed SN rates: ipsN based on the observed HII LF underestimates the observed value by an order of magnitude. The 
predicted value of ipsN leads to porosities < 1, similar to results for the other galaxies. However, the observed SN rate yields 
porosities S> 1, implying a strong dominance of the HIM. It is unclear how this discrepancy should be reconciled, since we 
have simply recast this pre-existing problem for the Galaxy. 

Our results suggest that OB superbubbles are indeed a dominant source of structure in the ISM of galaxies. We used the 
simplest and crudest formulations to derive the superbubble size distribution, with the aim of identifying dominant processes 
and evolutionary features. The agreement of this derived prediction with the high-quality SMC data is unexpectedly excellent 
and encouraging. This tentatively suggests that we have fairly successfully identified the basic features of shell evolution, 
ISM parameters, and stellar parameters that are relevant to the global structure and evolution of superbubbles. As new, 
high-resolution HI data become available for more galaxies (e.g., Thilker 1997; Kim etal. 1997), we will be able to more 
rigorously test these results and distinguish unresolved issues. 
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